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STABILIZATION OF TRANSVERSE VIBRATIONS OF AN 
INHOMOGENEOUS EULER-BERNOULLI BEAM WITH A 
THERMAL EFFECT 

OCTAVIO VERA, AMELIE RAMBAUD, AND ROBERTO E. ROZAS 


Abstract. We consider an inhomogeneous Euler-Bernoulli (EB) beam of length 
L clamped at both ends and subject to : an external frictional damping and 
a thermal effect. We prove the well-posedness of the model and analyze the 
behavior of the solution as t —>• -\-oo. The existence is proved using semigroup 
theory, and the exponential stabilization of solutions is obtained considering 
multiplier technique. A numerical illustration of the energy decay is given, 
based on initial data close to a real physical experiment. The results gener¬ 
alize the ones obtained where there is no heat conduction taken into account 
(see [a El [12] for EB type models, 11511161 for wave like equations), together 
with the coupled systems of type damped wave equation with heat conduction 

[SllTg]. 


Keywords: semigroup. Asymptotic behavior. Euler-Bernoulli. heat conduction, 
damping. Coupled system, multiplier. 

1. Introduction 

Simplified models of beams’ vibrations have interested researchers for decades. One 
of the main issues, both from a physical and mathematical point of view for such 
problems is the question of stabilization of the oscillations when time goes to inhn- 
ity. There exist various models of transverse vibrations of beams. Here, we focus on 
a Euler-Bernoulli (EB) like model, which is the most commonly used model because 
it is simple and provides reasonable engineering approximations. See for example 
[i m US [Ml El- However, the standard homogeneous EB equation can be gen¬ 
eralized in order to account for inhomogeneity of the material, heat conduction, 
and dampings (either viscous, or frictional). Let us also mention a model where 
both transverse and axial vibrations are investigated [26] • In the present work, we 
investigate a model of transverse vibrations, derived from the model proposed in 
[I]. Let L > 0 be the length of a axially symmetric thin beam, along the axis x. 
We also suppose that the microbeam is subjected to heat and to a non uniform 
frictional damping. Then the transverse displacement u and the temperature held 
9 satisfy the following set of equations (without dimension), for all t > 0, for all 
X G (0,L): 

( 1 . 1 ) utt + (pix) Uxx)xx + 2 q{x) Ut + n Oxx = 0 , 

( 1 - 2 ) 9i rjOxx i^Uxxt — 0 ; 

where p is a non negative function describing the non homogeneity of the beam, and 
q, also non negative stands for the force of the frictional damping. The constants 
K, rj are assumed to be strictly positive. Note that system (|n])-(IE2|) is a simplihed 
version of the one investigated in [1]. Indeed, although we added inhomogeneity 


1 


and damping, the authors in [I] consider a more general law of heat conduction and 
assume the beam is moving along the x axis at a constant velocity, and that the 
temperature is also damped. Moreover the coupling constants are not the same in 
their case (here we only have k > 0 ), so their problem is more difficult to tackle. 
System (frni-iLa is complemented by Dirichlet boundary conditions 

(1.3) u{0, t) = u{L, t) = 0, ^ 2 ,( 0 , t) = Ux{L, t) = 0, 0(0, t) = 9{L, t) = 0, 

that hold for all t > 0, describing a clamped beam at both ends. And finally we 
have initial conditions: 

(1.4) u{x, 0) = uq{x), Ut{x, 0) = ui{x), 6{x, 0) = 9o{x), x G [0, L], 

Our purpose in this work is to provide well-posedness of system (fO-lLll under 
smoothness assumption on the inhomogeneity functions p, q, and to investigate the 
long time behavior of the energy, that is to say wether or not the energy decays 
exponentially. Indeed, the study of the total energy functional £ : 1R+ —>■ 1R+, 
defined, in our case, for all t, by 

(1.5) £it) = i 

when t —>• oo, answers the question of the stabilization of the vibrations of the 
structure. Rates of decay for the total energy have been established for various 
dissipative models, see for example [312 El H il 13 [H [H] and references therein. 
Up to our knowledge, the exponential stability has not been investigated for our 
model, and provides a first answer to the stabilization of microbeams under study 
in [I], with of course only the Fourier law of heat, and a non moving structure in 
the X direction, for a start. Moreover, our result extends the results of USES] 
for inhomogeneous beams without heat, and the ones of |181119) for damped waves 
coupled with heat. 

The rest of the paper is organized as follows. In Section [3 we state the mathe¬ 
matical hypothesis and prove well-posedness of the initial boundary value problem 
(fOl-OTl) using a semigroup argument. In Section El we prove the exponential de¬ 
cay of the energy (11.51) which is the main result. Finally, in Section |2 we illustrate 
numerically the exponential decay of the energy by considering a simple Finite Dif¬ 
ference discretization, in the spirit of [ 22 ] and initial conditions close to a physically 
relevant situation, inspired by [23]. 

2. Well-posedness 

In this section, we prove existence and uniqueness of solution to (ini)-(iri|). Let us 
start by introducing the main assumptions and notations. 

2.1. Preliminaries. Throughout this paper, we always assume that q{x) and p{x) 
are positive definite functions, q{x) G L°°(0, L), p{x) G iL^(0, L) and that there 
exist constants ai, a 2 , and 0:4 such that 

0 < ai < p(x) < a 2 , VxG[0, L] 

^ ^ 0 < 03 < q{x) < 04 , Vx G [0, L] 

Throughout this paper (7 is a generic constant, not necessarily the same at each 
occasion(it will change from line to line), which depends in an increasing way on 
the indicated quantities. 
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We first establish an equality describing the dissipativity property of our system. 

Lemma 2.1. For every smooth enough solution of the system (HU-dLll) the total 
energy £ : M+ —>■ M+, defined by (11.511 satisfies 

d 

(2.2) —£{t) = -2 q{x)\ut\‘^dx-r] dx. 

Ctt Jo Jo 

Proof. We multiply (fTTl) . (HH) by ut and 9 respectively. Integrating in x G [0, L], 
using the boundary condition dOD, summing and taking the real part, the result 
follows straightforward. □ 

Remark 2.2. The equations (I1.5I) - (I2.2I) can be rewritten as 

^(0 = 2 (ll^/P^a:a;||i2(o, L) + ll^t ||i2(o, L) + II^1112(0. L)) 

and satisfies 

—£’(t) = — 2 ||^/9Mt|li2(Q_i) — r] ||0x|li2(o_L)- 

Let us now introduce the phase space associated to our problem. We use the 
standard L^(0, L) space, the scalar product and norm are denoted by 

pL pL 

(u, u)L 2 ( 0 _i) = UV dx, ||M||i 2(0 L)= dx. 

Jo Jo 

Hence, denote TL = Hq{0, L) x L^(0, L) x L^(0, L) endowed with the inner product 
given by 

pL pL pL 

{{u, V, 6), {ui, Vi, 9i))y^ = / p{x) UxxUixx dx + / vvidx+ / 9 9idx 

Jo Jo Jo 

for U = (tt, V, 9), U = (ui, vi, 9i) and the norm 

I|(u, V, 0)11-^ = \\'ypUxx\\L^(^o,L) + II'*^IIl2(o. L) + II ^*IIl2(o. L) • 

We can easily show that the norm || • ||^ is equivalent to the usual norm in "H, so 
that the phase space TL is a Hilbert space. 

Finally, we have the Poincare inequality 

ll^lli2(o. L) ^ C'p ||Mx||i2(o, L), VuG Hq{0, L), 
where Cp is the Poincare constant. 


2.2. Existence and uniqueness. In order to establish the well-posedness, we 
reformulate our problem as an abstract Cauchy problem. Taking ut = v, the initial 
boundary value problem (fO-dTHl) can be rewritten as: 

(2.3) fuit)=AUit), U{0) = Uo, Vt>0, 

dt 

where U{t) = {u, v, 9)'^ and Uq = {uq, Ui, 9fi)^. The linear operator 
A : ^{A) C "H —?> "H is given by 

V 

- 2 q{x) V - {p{x) Uxx)xx - K 9xx 
^ ^xx T rj 9xx 
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It remains to define the domain 'D{A) of the operator: 

V{A) = {(u, v,e)en: uG H^{0, L), - p(x) u,, - kO G H^(0, L), 

6»eij2(0, L), KV + 7]0Gff^(O, L)}. 

We are now ready to state our existence result. 

Theorem 2.3. For any Uq G H there exists a unique solution Uit) = (u, Ut, 0) of 
(or, equivalently, to the abstract Cauchy problem (12.3|) ). satisfying 

u G C'([0, oo[; Hllf), L)) (d oo[; ^^(O, L)) 

e G C([0, oo[; L'^iO, L)). 

However, if Uq G F>(A) then 

u G oo[; Hl{n)) n C2([0, oo[; L)) 

e G C'([0, oo[; Hl{0, L)) (d oo[; L)) 

P'^XX “1“ 1^0^ KUt + r]9 G C{[0, oo[; Hq{0, L)). 

Proof. Firstly, we show that the operator A generates a Co-semigroup of contrac¬ 
tions on the space H. We will show that is a dissipative operator and 0 belongs 
to the resolvent set of A, denoted by q{A). Then our conclusion will follow using 
the well known Lumer-Phillips theorem |20| . We compute: 

{AU,U).^= / p{x)VxxUa:x dx 
Jo 

rL 


+ / [-2 q{x) V - {p{x) Uxx)xx - K 6xx\ V dx-\ 

Jo 

pL pL 

2ilm / p{x) VxxUxx dx — 2 / q{x)\v\‘^ dx 

Jo Jo 


/' 


KVx -\-pOx)x 0 dx 


0 

nL p L 

2 inlm I Vx0x dx — p \9x\'^ dx. 

Jo Jo 

Taking real parts we obtain 

(2.5) Re{AU, C/)„ = - 2 ||y® - ry ||0.||i.(o,i). 

Thus ^ is a dissipative operator. Next, the domain 'D{A) is clearly dense in H 
and the operator is closed. Finally, let us show that 0 G p(A), the resolvent of 
A. In fact, given F = (/, g, h) G H, we must show that there exists a unique 
U = {u, V, 9) in ViA) such that AU = F, that is to say: 

(2.6) v = f in H^{0,L) 

(2.7) 2 qipf) V (p(x) Uxx)xx ^ 9 xx — 9 (65 C) 

( 2 . 8 ) KVxx + ridxx = h in L^( 0 , L). 

Taking v = f in (12. 8 |) we have 


(2.9) 


9xx = -(h-Kfxx) G L\0,L). 
V 


It is known that there is a nnique 9 G L) satisfying (|2.9|) and 


I1^xx||l2(0, L) ^ ^ 11^ ^fxx\\L^(0,L) ^ ^ 
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for a positive constant C. Moreover, we have that {kv + r] 6)xx = h G i^(0, L) and 
we conclude that {kv + r] 9) G Hq{0, L). 

Replacing (12.61) into (12.71) we have 

(2.10) {p{x)uxx)xx = - g -‘2q{x) f - kOxx in L^(0, L). 


Moreover 

rL 


/ dx= —tT p{x) ' 
Jo Jo 


da; < — 
ai . 


p{x) u 


dx= — \\V^ UxxWh^o^ 

CXi ^ 


L)- 


It follows using the Lax-Milgram theorem that there exists a unique function u G 
Hq{0, L) such that ()2.10p holds in L^(0, L). It is easy to show that WUWn < C ||F||h 
for a positive constant C. Therefore we conclude that 0 G p{A). Therefore A is 
the generator of a C° semigroup of contractions and the proof of Theorem 12.31 is 
complete. □ 


3. Asymptotic Behavior 


In this section we show that the energy decays uniformly exponentially with time. 
We state our main result as follows. 


Theorem 3.1. Let u, 9 he solution to system (fO-lfri). Then there exist positive 
constants M and m such that 

£{t) < M£{0)e-^*, Vt>0, 

where £ is the total energy of the system, given by CH). 


We will prove this theorem using multiplier technique and a suitable Lyapunov. 
We begin with preliminaries. 


3.1. Preliminaries. We recall the energy equality from Lemma [2T 


_d 

dt" jQ 

Let us now define an auxiliary function that will help building our Lyapunov: 


(3.1) 
Let u 

(3.2) 


pL/ nLi 

-£{t) = — 2 / q{x) dx — r] 9\ dx. 

1 Jo Jo 

auxiliary function that will help buildii 

p L/ pL 

Ai{t) = / uutdx+ q{x)u^dx. 


Jo Jo 

We conclude this subsection by a lemma to control the time derivatives of the 
auxiliary function Ai. 


Lemma 3.2. Let u,9 be a smooth solution to (fTB-lLl). Then the auxiliary func¬ 
tion defined by dH satisfies the following estimates. There exists p > 0 such 
that for all t > 0 : 

(3.3) — p £{t) < J-i{t) < p£(t). 

Moreover, the derivative of verifies: 

(3.4) = -2£{t)+ni{t). 

dt 

with the remainder TZi such that: there exist non negative constants Ci, C 2 , Co > 
0 , such that, for all a > 0, for any t > 0, it holds: 

dx+Cloa J pUxx- 


(3.5) 7^l(^) < Cl J + C 2 (^1 + J 
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Proof. We begin with the first estimate. From the definition of Ti, (13.21) . On 
the other hand, using the classical Young and Poincare inequalities successively, 
together with the assumption on functions p and q, dUD we obtain 

(3.6) < t: [ \ut\‘^ dx + C{Cp,ai,a 4 ) [ p{x)ul^dx, 

^ Jo Jo 

where C'(C'p,ai, 0 : 4 ) is a non negative constant, depending on Poincare constants, 
and the upper and lower bounds of p and p. The estimate (13.31) follows straightfor¬ 
ward with ^ 

/r = — C{Cp, ai, 0 : 4 ). 

Next, we differentiate Ti: by using (HH) and integrating by parts we have 


q{x)u dx 

pL 


*^ 1(0 — f uutt dx dx [ 

Jo Jo Jo 

= f u{- 2q{x)ut - {p{x)uxx) - f^0a;x)dxf uldx-\-^f q{x)u‘^ dx 
Jo Jo dt Jq 

pL p L pL 

= — p{x) dx — K 9 Uxx dx+ dx 

Jo Jo Jo 


— — 2£ -|- 1^1 , 
with 


7^l = 2 


Ut 


9^ -K 


9ux 


Hence, using the Poincare inequality, (12.11) and the generalized Young’s inequality, 
we have, for all a > 0 , for all t > 0 ; 

1^1 (01 < 2 / u^ + Cp f 9l + ^ f 9l + ^a ( p{x)ul^. 

Jo Jo 2 a Jq 2ai Jo 

Estimate (|3.5I) follows directly. We will chose a particular a > 0 in the following. □ 
With this Lemma, we are now ready to prove the exponential stability. 

3.2. Proof of Theorem 13.IL We introduce a Lyapunov functional: 

(3.7) g{t) = £{t) + e , 

where e > 0 is a positive constant to be adjusted later (as the a > 0 of the previous 
Lemma). 

Let us compute the derivative of g. From (13.11) . together with (13.41) and (13.51) from 
Lemma |T2l we get: 

§g{t) = -2j^^q{x)u^dx-rij^9ldx-2e£{t)+eni{t) 


(3.8) 


< - 2 - 


(2 - e Jo q{x) ui dx - (?7 - e C 2 (1 i)) 9l dx 

-£ (^2 £{t)- a C 3 Jq pu^x) ■ 


Now we will chose the constants e and a. We need to control the last term in (13.81) 
in order to have a good differential inequality. We also need to ensure the positivity 
of the Lyapunov g. 
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On the one hand, from the definition of the total energy, (CH), it holds: 

2 £{t)-C 3 a [ pul^ = (l-aCg) f pul^ + f u'i dx + [ \9\'^ dx. 
Jo Jo Jo Jo 

Hence, choosing a > 0 small enough so that 


K := l-aCo > 0, 

we get, injecting in (13.81) : 

dt 

where K > 0 is defined above, while a is now fixed to ensure K > 0, and: 

Ki = ^2 -e^^, K 2 = f] - eC2 (^1 + ^ . 

On the other hand, from the definition of G, dST]), and from (13.31) in Lemma 13.21 
we have, for alH > 0: 

(3.9) {1- e p) £{t) < G{t) < {1 + £ p) £{t). 

Therefore, it suffices now to choose e > 0, small enough so that: 

{ 1 — £ p > 0 , 

> 0 ’ 

K 2 = ?7 - eC2 (1 + i) > 0. 

Hence, the positivity of the Ki, i = 1,2 gives: 

(3.11) < -2K££{t), 

where £ is now fixed, so that (13.101) is satisfied. Finally combining (13.91) with (13.111) 
integrated in time, we get, for all f > 0: 


g{t) < - 2K££{t) - Ki 


q{x) dx — 


K2 


Oi dx 


£{t) < Mo- 


2K£ 


(1-£m) Jo 


£{s)ds, 


with 


Mn = 


^^( 0 ) 


l + £ p 


1 — £ p 1 — £ p 

The conclusion of Theorem [3JJ that is 


£:( 0 ) > 0 . 


£{t) < M£{Q)e-'^\ 

follows immediately from the Gronwall Lemma, with: 


2 K £ 1 + £ p 

m = z - >0, M = 


1 — £ p 


1 7 

— s ^ 


and the proof is complete. 
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4. Numerical simulation 


In this section we provide a numerical illustration of Theorem 13.11 that is the sta¬ 
bilization of the vibrations thanks to the thermal effect. Few numerical simulations 
have been don for a EB type equation, coupled with a damping and/or heat con¬ 
duction. See for example [22], where the authors use a Finite Element scheme in 
space, but do no care about the time discretization, since they are interested in 
the Optimal Control Model. Here we want to illustrate our result of exponential 
stability for physically relevant initial conditions. Therefore, since we do not aim 
at designing a high order scheme, we consider a classical second order Finite Dif¬ 
ference scheme in both time and space (see for example (Hj), taking care with the 
discretization of the bilaplacian operator ( [lain mi). We consider the following 
initial conditions: 


j u{x, 0) = 0, V a; £ (0, L), 

\ e{x,0) = Ae~^, V X £ (0,L), 

with A = 0.2, whereas the standard deviation changes from 10“^ to 10“^. 
These initial conditions represent a possible physical experiment, don to test the 
response of a beam (the deformation) when we apply an initial heat pulse in the 
center of the beam, see [23] . 

We plot in Figure |4| the discrete energy (computed with discrete integration at all 
discrete times of the continuous expression (11.51) 1. versus time in Log — Log 
scale, and observe the exponential decay for three values of the standard deviation. 

Conclusion 

In this work, we show the exponential stabilization of an EB type model of flexible 
inhomogeneous beam, taking into account frictional damping and heat conduction 
of Fourier type. We illustrate our result numerically, considering physically rele¬ 
vant initial conditions, and our result generalizes the results obtained for the wave 
equation and for EU models without heat. We think that this result is a first step 
toward a better understanding of more complicated models, such as the one derived 
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